Adiabatic nonlinear waves with trapped particles: I. General formalism 
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A Lagrangian formalism is developed for a general nondissipative quasiperiodic nonlinear wave 
with trapped particles in coUisionless plasma. The adiabatic time-averaged Lagrangian density £ 
is expressed in terms of the single-particle oscillation-center Hamiltonians; once those are found, 
the complete set of geometrical-optics equations is derived without referring to the Maxwell- Vlasov 
system. The number of trapped particles is assumed fixed; in particular, those may reside close to the 
bottom of the wave trapping potential, so they never become untrapped. Then their contributions 
to the wave momentum and the energy flux depend mainly on the trapped-particle density, as an 
independent parameter, and the phase velocity rather than on the wave amplitude a explicitly; 
hence, £ acquires a-independent terms. Also, the wave action is generally not conserved, because 
it can be exchanged with resonant oscillations of the trapped-particle density. The corresponding 
modification of the wave envelope equation is found explicitly, and the new action fiow velocity is 
derived. Applications of these results are left to the other two papers of the series, where specific 
problems are addressed pertaining to properties and dynamics of waves with trapped particles. 

PACS numbers: 52.35.-g, 52.35.Mw, 52.25.-b, 45.20.Jj 
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I. INTRODUCTION 

A standard approach to describing a nondissipative 
wave in the geometrical-optics (GO) limit is to start with 
its time-averaged Lagrangian, as proposed originally by 
Whitham in Ref. T; see also Refs. ^\^. Through that, 
both the nonlinear dispersion relation (NDR) and the ac- 
tion conservation theorem (ACT) are yielded, the latter 
being a particularly robust way to derive the envelope 
equation ITU]. However, the existing models using the 
time-averaged Lagrangian [51 [71 [TTlfT^] cannot account 
for effects caused by particles trapped in wave troughs. In 
particular, those effects require special treatment for they 
are not necessarily perturbative, i.e., may grow as the 
amplitude decreases [13]. Hence, describing waves such 
as Bernstein-Green-Kruskal (BGK) modes [TSHIS] has 
been limited to more complicated kinetic models [T71 [T^j- 
123] , which are specific to particular settings and may not 
render the underlying physics transparent. Therefore, it 
would be beneficial to generalize Lagrangian theories to 
accommodate trapped-particle effects. 

It is the purpose of this paper to do so. Specifically, 
a Lagrangian formalism is developed here for general 
nondissipative quasiperiodic nonlinear waves in collision- 
less plasma, under the assumption that the number of 
trapped particles remains fixed. This assumption is ob- 
viously satisfied for any stationary homogeneous wave. 
For nonstationary or inhomogeneous waves, the number 
of trapped particles being fixed implies that (i) those are 
trapped deeply, such that they do not become untrapped 
when the wave parameters evolve; (ii) also, it is implied 
that there are no passing particles close to the resonance, 
so that no additional trapping can result from the wave 
evolution. (As models, corresponding distributions al- 
ready proved useful for understanding paradigmatic ef- 
fects driven by trapped particles [24H27] : yet, they can 
also form naturally as waves evolve [10].) In particular, in 



the case of nonstationary or inhomogeneous waves, dis- 
tributions smooth across the resonance are not allowed. 

Under the aforementioned assumptions, we express the 
adiabatic time-averaged Lagrangian density £ in terms of 
the single-particle oscillation-center (OC) Hamiltonians; 
once those are found, the complete set of GO equations is 
derived without referring to the Maxwell- Vlasov system. 
Since the number of trapped particles is fixed within our 
model, their contributions to the wave momentum and 
the energy flux depend mainly on the trapped-particle 
density, as an independent parameter, and the phase ve- 
locity, rather than on the wave amplitude a explicitly; 
hence, £ acquires a-independent terms. (Of course, tak- 
ing the limit a —^ in £ would require that the width of 
the trapped-particle distribution also be zero, i.e., that 
the distribution be (5-shaped.) Also, the wave action is 
generally not conserved, because it can be exchanged 
with resonant waves of the trapped-particle density. The 
corresponding modification of the wave envelope equa- 
tion, or the ACT, is found explicitly for one-dimensional 
(ID) waves, a case in which the trapped-particle density 
is expressed directly in terms of the wave variables, thus 
providing an exact closure. 

The results presented here extend our Ref. [T3| in that 
we now (i) allow plasma parameters to vary slowly in 
space and time, and (ii) derive the corresponding en- 
velope equation, or the ACT, in addition to the NDR. 
Applications of these results are left to Refs. [551 US] 
(further referred to as Paper II and Paper HI), where 
specific problems are addressed pertaining to properties 
and dynamics of waves with trapped particles. 

The paper is organized as follows. In Sec.[ll| we derive 
the general form of £. In Sec. |III[ we consider ID waves in 
particular and obtain the corresponding ACT and NDR. 



In Sec. IV longitudinal electrostatic waves are studied as 
a special case. For an arbitrarily nonlinear wave, the par- 
ticle OC Hamiltonian is derived, generalizing the dipole 



ponderomotive Hamiltonian. Then, the action density, 
the action flux density, and the action flow velocity are 
inferred. In Sec. |V] we also calculate those quantities 
specifically in the small-amplitude limit. In Sec. |VI[ we 
summarize our main results. Some auxiliary calculations 
are also presented in appendixes. 



II. WAVE LAGRANGIAN 

In Ref. [14], we proposed the following expression for 
the Lagrangian spatial density of an adiabatic wave in 
coUisionless plasma |30j : 



(£c 



^n.(H. 



(1) 



Here (£om) is the time-averaged Lagrangian density of 
the electromagnetic field, summation is taken over dis- 
tinct species s, Us are the corresponding average densi- 
ties, and ("Hs)/^ are the corresponding oscillation-center 
(OC) Hamiltonians averaged over the distributions fs of 
canonical momenta p. The formula was originally de- 
rived for homogeneous stationary waves [31j . the case in 
which p and Ug are constants. What we show below is 
that Eq. (IT]) holds also in the general case, except now 
one needs to specify how Ug relate to the field variables. 



A. Plasma Lagrangian 



Consider the Lagrangian Ls — J £^ dV, with the 
spatial density £x; = ■Gem + ■Cp 32 . Here £cm = 
{E^ — B^)/{8n) is the field Lagrangian density, E ~ 
—Wifi — dtA/c is the electric field, B = V x A is the 
magnetic field, (p and A are the scalar and vector poten- 
tials, and c is the speed of light. Also, 



£p = } ^(5(x-Xj)£i(x,v,;y), A), 



(2) 



where the summation is taken over individual particles, 
and Li are the Lagrangians of those particles; namely, 
Li = LI ' + L^ , where L| are independent of the 
field, and L^ ' — (ei/c)(vj • A) — Ciip. Finally, Xi(i) are 
the trajectories of individual particles, v^ = x^ are the 
corresponding velocities, and e^ are the particle charges. 
Suppose that the electromagnetic field contains a 
rapidly oscillating part and consider the plasma dynam- 
ics on scales large compared to the oscillation scales. (In 
the presence of resonant or trapped particles, one of such 
scales is the period of bounce oscillations Tf, [33 Sec. 8- 
6].) Then, it is only the time-averaged part of the La- 
grangian, Ls = /(£s) dV, that contributes to the system 
action. Hence Ls plays a role of the slow-motion La- 
grangian of the system pT . Specifically, we write 



(iiE> = (iiom) + (4'')) + (4*'). 



Here (£om) generally consists of two terms, £cin due to 
quasistatic fields {(p,A) (if any) and (£em) due to the 
actual wave field. The remaining terms describe contri- 
butions of passing particles and trapped particles, corre- 
spondingly, and are derived as follows. 

In the case of passing particles, we separate the slow, 
OC motion Xi(t) and the quiver motion Xi(t) and notice 
that 

y(5(x-x,-x0i[''^(x,v,))dV = 

y"(J(x-x,)£l^^(x + x„v,-fv,))rfV = 

j 5{^ - i,) (l|p^(x + i,, V, + v,)) dM = 

y'<5(x-x,)4'')(x,v,)dV. (4) 



Here we introduced 

4p^(x,v) = (L,(x + x„v + v,)), 



(5) 



which has the meaning of a single-particle OC Lagrangian 
[34) . Hence, one can write 



(4^))^^<5(x-x,)£|^)(x,v,) 



(6) 



In the case of trapped particles we proceed similarly, 
except that the OC location Xi(i) is now determined by 
the motion of the wave nodes and thus cannot serve as an 
independent variable. Instead, the new independent vari- 
able will be the phase Oi of bounce oscillations, possibly 
in multiple dimensions. Since these bounce oscillations 
are assumed adiabatic, Cl will not depend on Oi explic- 
itly; rather it will depend on Oi and, parametrically, on 
Xi. (Remember that the dependence on the field vari- 
ables is also implied throughout the paper.) Thus, the 
Lagrangian of the bounce motion, henceforth also called 
OC Lagrangian for brevity, is given by 



£f^(x,0,) = (L,(x + x„u + v,)), 



(7) 



where we substituted the wave phase velocity u for the 
average velocity. This yields 



(4)) = ^%-x,)/:f(x,0O 



(8) 



B. Routhian 

Below, it will be more convenient to use canonical OC 
variables for particles, (q^, Pj). For simplicity, let us tem- 
porarily require that q^ = x^ for passing particles and 
qj = Oi for trapped particles; in the latter case the canon- 
ical momentum Pj will be the action J^ of the bounce 
oscillations. Then, let us use [35] 



(3) 



C. 



PiQi 



H. 



(9) 



so Eq. (pl rewrites as 



where L = / £ dV, and 



L^ 



J(x-x,)-Hi(x,p, 



(10) 



(11) 



Since both q, (t) and Pj (i) are now independent functions 
(cf. Ref. [36l Sec. 40]), field equations wih be insensitive 
to the second term in Eq. (10 1; i.e., for the purpose of 



finding jield equations, this term can be dropped. There- 
fore, L plays the role of the adiabatic Lagrangian of the 
wave (and also of quasistatic fields, if any). 

Notice that, since L = Lj; — X^i Pi ' ^n i* can be con- 
sidered as a Routhian of the particle-field system [5S1 
Sec. 41], i.e., a function that acts as a Lagrangian for 
the field variables but as a Hamiltonian for the particle 
variables. As a Routhian, the wave Lagrangian was also 
introduced earlier in our Ref. |14j . (However, unlike in 
Ref. |14|, here we do not perform Routh reduction (Ap- 
pendix Kl); i.e., now we allow p to evolve.) Below, we 
will show how Eq. ( 11 ) corresponds to that earlier result. 



C. Locally averaged densities 



In Eq. (11), the summation over all particles i can be 



separated into (i) summation over species s, (ii) sum- 
mation over p ■ within a local elementary spatial volume 
AVfc, and (iii) summation over all AV^. Specifically, let 
us choose the elementary volumes large enough such that 
both Hs and the densities n^ vary little fST" within AV^:. 
Then, 

^5(x-x,)Hj(x,pJ 

i 

^^^^{^-^k) ^■H,(xfc,Pj) 

s k j 

= X! X! ^(^ ^ ^^^ ris(xfc, ^)('Hs(xfc, p))/^ AVfe 

s k 

= ^n,(x,i)(H.(x,p));^. (12) 

S 

This puts £ in the anticipated form, Eq. ([l]), or 

s s 

From now on, the specific canonical variables will not 
matter; i.e., further canonical transformations are al- 
lowed in Hs and "Hs , if necessary. 

D. Independent variables 

We are now to choose the independent variables that 
will describe the field. [Quasistatic fields, if any, can be 



described by (i^. A) as usual and thus will not be con- 
sidered explicitly] Suppose that the wave is character- 
ized by a smooth envelope a(x, i), arbitrarily normalized. 
Also suppose that the wave field, while not necessarily 
monochromatic, oscillates rapidly with some canonical 
phase ^, the period being 27r. Hence, the local tem- 
poral and spatial periods can be defined as T = Itx jio 
and A = 2t: /k, such that u = uj/k is the phase speed, 
u = ?ik/fc. 



io = -dt^, k = vc, 



and, in particular. 



dth + diio = 0, 



Ojki 



dikn = 0. 



(14) 



(15) 



where we introduced di = d^^ ■ The function £ will then 
depend on {dtC ^0 E^. b^t ^^t on ^, for it describes 
the dynamics on scales t and A such that t ^ T and 
£ >• A Ij. (In the presence of trapped particles, we also 
require t ^ ti, and £ ^ utj,.) 

The question that remains is how to treat n^ when 
varying £. For passing particles, the OC densities nl 
are determined by x, = q^, which are independent vari- 
ables; thus, nj' (x, t) are also independent of the field 
variables. For trapped particles, however, Us (x, t) are 
determined by x^ which may be tied to the wave troughs; 
thus, ris (x, t) may be connected with the wave phase. 
In particular, for ID waves this connection can be imple- 
mented as an exact closure, which is done as follows. 



III. ONE-DIMENSIONAL WAVES WITH 
TRAPPED PARTICLES 

A. Extended Lagrangian 

First of all, notice that, in a ID system, trapped par- 
ticles travel at the wave phase velocity u. Hence, the 
corresponding continuity equations read as 



dtns +dx{nsu) = 0, 



(16) 



where we introduced n, = nL to shorten the notation. 



One can embed Eq. ( 16 1 in the formalism by considering 

(17) 



a new, extended Lagrangian density 

A = £ + ^ ^s [dtns + dxixisu)] . 



Here /ig are Lagrange multipliers |39| . i.e., new indepen- 
dent functions of (i,a;), yet to be found. In particular, 
varying A with respect to ris yields 



dttxs + ud^^is + {n<-/^)f^^o, 



(18) 



whereas Eq. ( 16 ) flows from varying A with respect to /i^. 

(19) 



Further, notice that 

A = £ - ^ ris [dtfis + i^/k) d^fJ-s 



is a Lagrangian density equivalent to A, yet with an 
advantage that A depends on only [38] the first-order 
derivatives of ^ [cf. Eq. (14)]: 



A = A(a, dtt d^S,, n, St/z, d^n). 



(20) 



Then, varying A with respect to ^ is as usual and yields 
dtkuj — dx^k = [T]. (We henceforth use indexes uj, k, 
and a to denote the corresponding partial derivatives.) 
On the other hand, 



A^ = 2u, - T^^s dxf^s, 



Afe = £fc + — y xis dxfis, 



(21) 
(22) 



where the derivatives are taken, in particular, at fixed n. 
Thus, one obtains 

9t£^ - a,£fc - ^ M„ (23) 

S 

Ms = dt{(Ts dxUs) + dx{(Tsud.j;fJ.s)- (24) 

Here we introduced as = x\s/k, which is proportional to 
the number of trapped particles (of type s) within one 
wavelength. For adiabatic waves this number is constant 
in the frame moving with the phase velocity; i.e.. 



dtCTs +udx(Ts = 0, 



(25) 



which is also obtained from Eqs. ( [T5| and (16). (No- 
tice that the equation for Us does not have a form of a 
continuity equation, unlike that for ris.) Then, 

Ms = as [^l^^is + ^x{u^.,^ls)\ - -^s 9.(^1*))/,, (26) 
where we used Eq. ( jlS] ). Hence, one gets 



B. Action conservation and wave dispersion 



From now on, let us consider £ as [38] 
£ = £(a,9t^,a,e,^). 



Using Eq. (29), one can hence write Eq. (28) as 



dtZ^u, — dxS,k 



/^ 



{ni'^)f.dx<ys 



(30) 



(31) 



Equation (31) represents a generalization of the well- 



known ACT for ID waves fTj, reproduced in the limit 
as — (also see Appendix^. Thus, we interpret 



2: = £. 



J = -Qk 



(32) 



as the new action density and the new action flux density. 



correspondingly. Notice, however, that Eq. (31), or 



dtI + dxJ^Y.(H 



(t) 



)fs dxC^s 



(33) 



does not have a conservative form in general, due to the 
nonzero right-hand side. This is because the wave of the 
trapped-particle density is, by definition, always resonant 
with the electric field, so the two can exchange quanta 
whenever as are modulated. Interestingly, the effect of 
cr- waves, which are described by Eq. ( |25[ ), is similar to 
the effect of entropy waves on magnetohydrodynamic os- 
cillations reported in Refs. [iOlUT] . 

Finally, we can complement the ACT with the NDR, 
by varying £ with respect to the wave amplitude. With- 
ing the lowest-order (in t~^ and ll~^) GO approximation, 
to which we adhere throughout the paper, £ depends on 
the local amplitude a but not on its derivatives [Eq. (30)], 



as usual [T, "47 . Remember also that as are introduced 
as independent functions and thus do not depend on a 
by definition. Hence, the NDR attains the same general 
form as for adiabatic waves without trapped particles, 
namely. 



£a=0. 



(34) 



or, equivalently. 



ft£. 



dx 



-£k + J2<Js{n 



(t)\ 



Y.{ni'^)f^dxas. (28) 



Further notice that {£,k)a = (£fe)n + J2k crs{dn,&)k, so 
(£.). = (£.)n-E^«(^«'^)/- (29) 



IV. LONGITUDINAL WAVES 

Now let us consider the specific case of ID longitudinal 
waves as a paradigmatic example. To do so, we will need 
to construct the Lagrangian density £ for such waves, 
which, in turn, requires calculating the single-particle OC 
Hamiltonians Us first. For linear waves, this is done 
in Appendix |B] For nonlinear waves, this also can be 
done straightforwardly, at least to the zeroth order in 



and 



Namely, we proceed as follows. 



where we substituted n^ = ask, in the left-hand side; 
similarly, {2,uj)<t — (£t^)n- (Here the external subindexes 
show variables kept fixed at differentiation.) Thus, it is 
convenient to consider £ as a function of as rather than 
of n^, specifically as follows. 



A. Single-particle OC Hamiltonians 

Consider the Lagrangian of a single particle in a sta- 
tionary homogeneous electrostatic wave in nonmagne- 
tized plasma. It is only the ID motion along the wave 



field E that matters for us; thus, we take 
L = mv /2 — e(p{x — ut), 



(35) 



where m and e are the particle mass and charge, v is the 
velocity in the laboratory frame K, and the potential (p 
is periodic yet not necessarily sinusoidal. (A quasistatic 
potential (f can be included straightforwardly and will 
not be discussed here explicitly.) Rewrite L as 

L = mu /2 + muw + mw" — £, (36) 

where y — x ~ ut, so y = w ~ v — u, and also 

,,2 , 



£ — mw /2 + e(f(y), 



(37) 



which is the energy in the moving frame K where (f is 
static. Then, time-averaging yields 



(L) = mu /2 + mu{w) + m{w ) — £, 



(38) 



where we used that £ is conserved on the oscillation scale. 

To proceed, it is convenient to introduce the angle 9 

and the action J of the oscillations in K; in partucular. 



J — -— (b w du. 

27r /^ ^ 



Then the oscillation period can be expressed in terms of 
the corresponding canonical frequency f2, yielding {w) = 
q^jk, where <^ = sgn (w) — sgnw (for trapped particles 
<; = 0), and m(u;2) = jn. Then, 



(L) == muV -mu^/2 + £, 



(40) 



where we used that the average velocity V = (v) equals 

V :^u + qQ/k (41) 

and introduced £ = JQ — £ . 

In particular, notice the following. Since the gener- 
ating function of the transformation (y, mw) — > (0, J) 
clearly does not depend on time explicitly, £ acts as a 
Hamiltonian in (0, J)-representation (cf. Ref. [H]), and 
thus £ is the corresponding Lagrangian; then. 



dj£ ^ n, dj£ ^ Jdjn 



(42) 



[Partial derivatives are used because £ and £ can also 
depend parametrically on {a,uj,k); cf. Sec. IV B[ ] 

First, let us consider a passing particle. In this case, 
for the canonical momentum p one can take P = dv{L) 
[Ml H5] . Assuming J = J{V), one thereby obtains P — 
mu + dj£dvJ, or, using Eq. (42 1, P = mu + J9ySl. 
Hence, from Eq. (41), we get 



P = mu + qkJ. (43) 

Then the OC Hamiltonian Ti. = PV — (L) reads as 

n^£ + Pu-mu^/2. (44) 



In particular, notice that Eq. ( 44 1 can be understood as 



a generalization of the nonrelativistic dipole ponderomo- 
tive Hamiltonian [Eq. ( [bT| ), with $5 from Eq. ( |B13| )] to 
the case of fully nonlinear particle motion in an arbitrary 
longitudinal electrostatic wave. 

In case of a trapped particle, the average coordinate 
is fixed, yielding T^ = u, so now those are {9, J) that we 
choose to serve as (q, p). Hence, V. — Jfl— (L), and thus 



'H^£ - mu'^/2. 



B. Parametrization. Wave Lagrangian 



(45) 



Although a can be defined as an arbitrary measure of 
the field amplitude, for the purpose of this paper it is 
convenient to introduce it specifically as the amplit ude 
of the wave electric field E. Hence, we can write Eq. ( 13 ) 
explicitly as 



167r 



-j:^i^Hni^^)j^-J2ni'Hn 



(t)N 



fs 



(46) 



with Hs to be taken from Eq. (44), and Us to be taken 



(39) fromEq. (45) 



Notice also that, once we have adopted a = E, the 
bounce-motion energy £ can depend parametrically on a 
and k but not on w, because the particle motion in K is 
entirely determined by the spatial structure of the wave 
potential, which is static there. In other words. 



£ — £{J,a, k). 



(47) 



Yet note that J — J{P,uj,k) for passing particles [see 
Eq. (43l], whereas for trapped particles J is an inde- 



pendent variable. In particular, this yields the following 
equalities that we will use below. First of all. 



d^£<-P^ = dj£'-P^ 5„ J(P) = -m<,Vtlk^ 



(48) 



where we substituted Eq. ( 42 ) for dj£ and Eq. ( 43 1 for 
di^J{P,uj,k); also, di^£'^*'^ — 0. Hence, one obtains 

du.n'^P'' = {P- mV)lk, a^-H^*' = -mu/k. (49) 



C. Action density 

Now we can calculate the wave action density T 
[Eq. ( |32|], namely, as follows. Since the first two terms 
in Eq. ( |46[ ) are independent of lo and k, one gets 

2: = Z(P)+Z(*), Z('')=-^ni')(a^Hf)/,, (50) 



(51) 
(52) 



with b = p,t. Hence, Eq. ( |49[ ) yields 



l(*)=fc-i^ 



rv-^'mM. 



One may say that X^^^ is proportional to (the density 
of) the ponderomoti ve rn omentum carried by passing par- 
ticles (cf. Appendix B 2 ), and I'-*-' is proportional to the 
kinetic momentum carried by trapped particles. Also, 
notice that the two can be combined as 



fcZ = ^n,m3(y)/^-^n(P'(P) 



/=• 



(53) 



The right-hand side here equals the difference between 
the system total kinetic momentum [ug = ns + ns 
being the total density of species s] less the momentum 
stored in particles, i.e., that of the untrapped population. 
Therefore, by definition, kX represents the wave total 
momentum, in agreement with Ref. [31 Sec. 15.4] (see 
also Paper III). Notice that a part of this momentum 
[namely, fcl*^*^] is independent of a, because it is stored 
in the trapped-particle translational motion with velocity 
u — Ld/k. However, remember that we still assume that 
a must remain large enough, such that detrapping does 
not occur (Sec.|l]). 

D. Action flow 



The action flux density J' [Eq. (32)] can be found sim- 
ilarly from Eq. (46) and reads as J' = J'^p^ + J^^\ where 



^W=^nWa,(H 



(b)^ 



/=' 



(54) 



again with h = p,t. In particular, notice that a part 
of J is also independent of a, like X^*) . This is because 
J = n/w, where 11 is the energy flux density |H Sec. 15.4] 
(see also Paper III) , a part of which is determined by the 
trapped-particle average velocity u rather than a. Notice 
also that vx = J^/I has the meaning of the action flow 
velocity and reads as 



In particular, with g 
one gets 



lip) +i(t) ' 



(55) 



Z(*)/I(rt and 4"' = J('''i/I^'\ 



,(p) 



vi = 



+ QVi 



(0 



1 + Q 



(56) 



V. DISCUSSION 

Now that we have developed the general formalism, 
it is instructive to consider ID longitudinal electrostatic 
wavesat small a in particular. Then, one can use 
Eq. (|bi| for v}f^ 



£(0) 



, SO 

167r 



Y.^'^\n^')L. 



(57) 



where U'-"^ is independent of the field variables, and e is 
the longitudinal dielectric function (Appendix IbI). Also, 
l(p) w ei^\E\^ / {Xdm) , and thus wf* equals the linear 
group velocity VgQ — ujk (Appendix B2|. Further, let 



At) 



us neglect £s compared to rrisU /2 in Eq. (45). Then, 



Hs ~ — msW^/(2fc^), so one obtains 



(58) 



which is independent of a, as expected. This yields 
vj- ~ u, and, therefore. 



Vx 



VgQ + 0U 

l + g 



(59) 



Equation ( 59 ) should not be confused with a similar 
expression in Ref. [15] derived for what is called there 
the nonlinear group velocity. The effects addressed in 
Ref. [45] result in a nonconservative form of the envelope 
equation, i.e., violate the ACT; hence, they are dissipa- 
tive by definition. In contrast, our formulation does not 
account for collisionless dissipation (except at inhonio- 
geneous as); thus, in Eq. (59) the difference between vx 
and Vgo is entirely due to adiabatic effects [¥7] . 

In addition to vi that we presented, one can also define 
the velocities of the energy and momentum flows [T] . For 
a nonlinear wave, those will be different from each other 
and from the true nonlinear group velocities Vg, at which 
the modulation impressed on a wave propagates adiabat- 
ically [15] . The effect of trapped particles on those true 
Vg will be discussed in Paper III. 



VI. SUMMARY 

In this paper, a Lagrangian formalism is developed for 
general nondissipative quasiperiodic nonlinear waves in 
collisionless plasma. Specifically, the time-averaged adia- 
batic Lagrangian density is derived in the following form: 



£=(iiom)-5I"-(^^)/= 



(60) 



Here (£om) is the time-averaged Lagrangian density of 
the electromagnetic field, summation is taken over dis- 
tinct species s, Ug are the corresponding average densi- 
ties, and {T-Ls)/, are the corresponding OC Hamiltonians 
averaged over the distributions fs of canonical momenta. 
Once Hs are found, the complete set of GO equations is 
derived without referring to the Maxwell- Vlasov system. 
For the first time, the average Lagrangian accounts also 
for particles trapped by the wave, under the assumption 
that the number of these particles remains fixed (Sec.ll|. 
In particular, ID waves are considered, in which case £ = 
£(a,a;, fc, cr); here w and k are the wave local frequency 

and the wave number, and ag = ris /k are proportional 
to the number of trapped particles within a wavelength. 



Correspondingly, the GO equations are summarized as 
follows. The first one is the consistency condition, dtk + 
dxtjJ = 0, due to cj = — 9t^ and k — d^S,; here ^ is the 
wave canonical phase. The second one is the NDR, given 
by £a = 0. The third GO equation is a modified ACT, 



dtl + d,J 



E 



W'pd^as, 



(61) 



with Z = £tj and J^ = — £fc being the action density and 
the action flux density, correspondingly. Because of the 
source term on the right-hand side, the wave action may 
not be conserved then, due to the fact that it can be 
exchanged with resonant waves of the trapped-particle 
density (u- waves). 

Since the number of trapped particles is fixed within 
our model, their contributions to the momentum density 
p and the energy flux density 11 depend mainly on Us and 
the phase velocity u, rather than on the wave amplitude a 
explicitly. Hence, T and J7 (and thus £ too) may contain 
a-independent terms, because 



I = p/k, J = n/w. 



(62) 



Particularly, in the limit of small a, the action flow ve- 
locity is obtained. 



«i 



1- 



- QU 

Q 



(63) 



where Vgi^ = Wfe is the linear group velocity, and q = 
X(*)/X(p)_ The difference between vi and w^q here is en- 
tirely due to adiabatic effects, so Eq. (63) should not be 



confused with a seemingly akin formula in Ref. j45| . 

Applications of these results are left to Papers II and 
III, where specific problems are addressed pertaining to 
properties and dynamics of waves with trapped particles. 
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Appendix A: Routh reduction 

In this appendix, we restate the concept of Routh re- 
duction |48) . complementing the derivation of the wave 
Lagrangian that we reported earlier in Ref. [13] . 

Consider a dynamical system described by generalized 
coordinates q = {q^,. . .q^), so the corresponding La- 
grangian has a form L(q, q, t), t being the time. The 
original least-action principle is then formulated as fol- 
lows 3.6] Sec. 2]: among trajectories q(i) starting at q^ 
at time ti and ending at q2 at time i2i realized is the one 
on which the action S = J^^ L dt is minimal. Since the 



general variation of S reads as 



(55 = 



dL 
dq 



d dL 
dt 9q 



Sqdt + p ■ Sq 



ti 



(Al) 



with p = dqL, and (5q(ti_2) = due to q^ 2 being fixed, 
one thereby obtains the Euler-Lagrange equations 



d^dL 
dt dq 



dL 
dq 



(A2) 



For simplicity, we henceforth consider a 2D system, 
with q — {9,x), the extension to a larger number of di- 
mensions being straightforward. Suppose, in particular, 
that 6 is cyclic, i.e., does not enter L explicitly. Then the 
momentum J canonically conjugate to 9 is conserved, so 
we can use the equality J = dnL to express as 



r2(J, X, x,t). 



(A3) 



Hence, S can be understood as a functional of x{t) only. 
Now consider the set of trajectories x{t) starting at 
given xi at time ii and ending at given X2 at time ^2, 
while 01 and 02 are arbitrary [albeit connected through 
Eq. (|A3|)]. Suppose that x{t) satisfies Eq. (A2| and 



consider the linear variation of S[x(t)] with respect to 
Sx{t) = x{t) — x{t). Then, from Eq. (|A1[), one obtains 



6S^Ji 



(A4) 



Since S is thereby not minimized on x{t) within these 
variation procedure, consider another, "reduced" action 



S = S 



J0dt. 



(A5) 



Here, the latter term can also be put as /^^ J dO^ where 
depending on x{t) are only the integration limits. There- 
fore, its variation around x{t) equals J 50\\\, thus yield- 
ing ^5'[a;(i)] = 0. Then a new variational principle can be 
formulated as follows: among trajectories x{t) starting at 
xi at ii and ending at x^ at t2, with arbitrary 0i and 02, 
realized is x{t) on which 5* is m inimal. 

Notice further that Eq. (A5) rewrites as 5 = J^ ^ Rdt, 
with the equivalent Lagrangian 



R{x, X, t) = L{x, X, t, n) - jn 



(A6) 



[where 51 — fl{J,x,x,t)], also known as Routhian. Hence, 
the motion equation that flows from the new variational 
principle reads as 



d dR dR 
dt dx dx ' 



(A7) 



also in agreement with the general Routh equations |361 
Sec. 41]. Since R is independent of 9, the system phase 
space is effectively reduced now, and a;-motion decou- 
ples, which is what constitutes the Routh reduction. 
For plasma physics applications of this technique, see 
Refs. Em 



Appendix B: Linear waves 



The quantity I = ii^j can be written as 



Here, we show how the known GO equations for hnear 
electromagnetic waves [which have n'-*^ — 0] follow from 
the general Lagrangian formalism discussed in Sec. ITT] 



1. General electromagnetic waves 

Let us take E, B ex e~*'^*+'''''', and employ the dipole 
approximation for H^ [49l [50] , namely, 



Us = nf^ + $s, $. = -E* • 6c, ■ E/4, 



,(0) 



(Bl) 



where "Hs is some function of the particle canonical mo- 
menta, $s is the ponderomotive potential, and ccg is the 
linear polarizability. Then, since 



l + 47r^ns(Q:s)/^ = e(w,k). 



where e is the linear dielectric tensor, one obtains 



where the term 



IGtt 



E* 



•E-IBI 



Y,ns{n 



(0)\ 



(B2) 



(B3) 



(B4) 



is independent of the wave variables. Now let us intro- 
duce the wave amplitude a via E = ae for the electric 
field envelope, where e determines polarization; hence, 
_B = |n X e|a, where n = ck/w, and c is the speed of 
light. Then, 



£ - £(") 



IGtt 



S(w,k), 



where we introduced 

2)(w,k) 



(B5) 



(B6) 



Varying the wave Lagrangian with respect to the am- 
plitude a yields the dispersion relation £„ = 0, or 



S)(w,k) ==0, 



(B7) 



which coincides with the known dispersion relation at 
prescribed e [33l Sec. 1-3]. In fact, the vector equation. 



e • e + n X (n x e) = 0, 



(B8) 



also can be recovered, namely, by varying £ with respect 
to e*. [One coul d, o f course, vary £ also with respect to 
E* and get Eq. ( B8 ) immediately.] 

Now let us vary £ with respect to the wave phase ^, 
with Eq. ( )l4| taken into account. Like for any other 
Lagrangian density of the form £(a,a;,k), one obtains 
then \^ 



dt£u. - V • £k = 0. 



(B9) 



167r ^ 



IGttlj 
1 

IOtto; 



3* • (cje^) • e + e* • e • e 
E*-9^(eL^)-E-f |B|2 



(BIO) 



or I = e/w, where we used Eq. (B7) and introduced e 
for the linear- wave energy density [49 j thus, I equals the 
linear-wave action density. Also, JT" = — £k can be put as 



J 



167r 



®i 



167r 



2). 






VgO^r, 



(Bll) 



where we used —d\^'D / du^l) — cjk [from Eq. (B7)], the 



latter being the linear group velocity Vgo; thus, J' is 
the linear- wave action flux density. Hence Eq. (B9) 
rewrites as 



dtl + W- (v„oI) - 0, 



(B12) 



in agreement with the linear ACT [2, Sec. 11.7]. 

Equations (14|, (B7), and (B12) represent a complete 



set of equations describing nondissipative linear electro- 
magnetic waves in the GO approximation (cf. Ref. |U 
Chaps. 14, 15]). As one can see from the above calcu- 
lation, the Maxwell's equations and the Vlasov equation 
per se are not needed to derive these equations 53]. 



2. Longitudinal electrostatic w^aves 

Finally, let us consider longitudinal electrostatic waves 
in somewhat more detail. In this case, for the longi- 
tudinal polarizability of an individual particle with OG 
velocity V, we take a^ = — e^/[ms(w — k- V)^] [33], where 
Cs and nis are the particle charge and mass, respectively; 
in particular, this corresponds to 



$., = 



e lEj^ 



(B13) 



(cf. Refs. 
function. 



4ms (w - k • V)2 
IIMIISS]). Then, the longitudinal dielectric 



e = 1 + 47r^7is(as)/. 



can be written as (cf. Ref. [53]) 

/s(T4) 



e=l-Y,^psj_ 



-oo (^ - kv^y 



dV^, 



(B14) 



(B15) 



where uj^^ — Anrise'^/ms, and fs{Vx) are the distributions 
of the particle longitudinal velocities Vx, normalized such 
that /_ fsiVx) dVx — 1. By definition, a linear wave has 
no trapped particles, so fsiVx) is zero in the resonance 



vicinity, and thus the integrand in Eq. ( B15 1 is analytic. so X = k ^ J^s "s i^V-^ *^s) /s • From Eq. ( B13 1 , one gets 



Hence, one can take the integral by parts. This yields 



= '+T.^-v 



LJ — kVr 



dV,, 



(B16) 



in agreement with Ref. [56 . 

From, Eq. (B7| the dispersion relation now reads as 
e(a;,k) = 0. In particular, this means v^o = — Ck/ecj and 



1 = 



IGtt 



lEl' 



(B17) 



Let us show that this expression is consistent with 

(B18) 



Eq. ([51]). First, notice that f49| 



fc-i(%$ 



S}fs 



^{d^^s)f^^-d^{^s)f^, (B19) 



and therefore I — —d^j X^s ''^s{^s)f^, or 



|-g|2 



(B20) 



where we substituted Eq. (Bl) for $3. Using Eq. (B14|, 
one thereby matches Eq. (B17), as anticipated. 
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